Confinement and Superfluidity in one-dimensional Degenerate Fermionic Cold Atoms 
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The physical properties of arbitrary half-integer spins F = N — 1/2 fermionic cold atoms trapped 
in a one-dimensional optical lattice are investigated by means of a low-energy approach. Two 
different superfluid phases are found for F > 3/2 depending on whether a discrete symmetry is 
spontaneously broken or not: an unconfined BCS pairing phase and a confined molecular superfluid 
instability made of 2N fermions. We propose an experimental distinction between these phases for 
a gas trapped in an annular geometry. The confined-unconfined transition is shown to belong to the 
Zn generalized Ising universality class. We discuss on the possible Mott phases at 1 /2N filling. 
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In the recent past, spectacular experimental progress 
has allowed to cool down alkali atoms below Fermi tem- 
peratures offering the promising perspective to study 
strongly correlated electronic effects, such as high- 
temperature superconductivity, in a new context 0. On 
top of strong correlations, ultracold atomic systems pro- 
vide also an opportunity to investigate the effect of spin 
degeneracy since in general alkali fermionic atoms have 
(hypcrfmc) spins F > 1/2 in their lowest hyperfinc man- 
ifold. Simplest examples are spin-3/2 atoms such as 
9 Be, 132 Cs, and 135 Ba atoms. In this Letter, we eluci- 
date the generic physical features of half-integer spins 
F = N — 1/2 fermionic cold atoms in the particular 
case of a one-dimensional optical lattice. The low-energy 
physical properties of these systems are known to be de- 
scribed by a Hubbard-like Hamiltonian Q : 



H = -t 
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where c s ai (a = 1,...,2N) is the fermion creation oper- 
ator corresponding to the 2F + 1 = 2N atomic states. 
The pairing operators in Eq. Q arc defined through 
the Clebsch-Gordan coefficient for two indistinguishable 
particles: Pj Ma = ^ ap (JM\F, F; a,f3)cl/ % . The in- 
teractions are SU(2) spin-conserving and depend on Uj 
parameters corresponding to the total spin of two spin- 
F particles which takes only even integers value due to 
Pauli's principle: J = 0, 2, 2JV — 2. Even in this 
simple scheme the interaction pattern is still involved. 
It is thus highly desirable to focus on a much sim- 
pler paradigmatic model that incorporates the relevant 
physics of higher-spin degeneracy. In this respect, we 



consider a two coupling-constant version of model l(T]l 
with U 2 = ••• = U 2 n-2 7^ U a : 

W = -*Z)[4 I< Ca,m+H.C.] + T X)^ +y S 7r i 7ri ' 



Pi = c a jCa.,% being the particle number operator and 
U = 2U 2 , V = U - U 2 . In Eq. ©, the singlet BCS 
pairing operator for spin-i* 1 fermions is 7r| = \/2N Pq i = 

c a iJapCg it the matrix J being the natural generaliza- 
tion of the familiar antisymmetric tensor e = io~ 2 to spin 
F > 1/2. Such a singlet-pairing operator has been ex- 
tensively studied in the context of two-dimensional frus- 
trated quantum magnets Q. When V = 0, i.e. Uq = U 2 , 
model (J2J) corresponds to the SU(2iV) Hubbard model. 
The Hamiltonian @ for V / still displays a large sym- 
metry since it is invariant under the Sp(2iV) group which 
consists of unitary matrices U that satisfy U* jW = J 
4J. In the spin F = 1/2 case, model © reduces to 
the SU(2) Hubbard chain since SU(2) ~ Sp(2). In the 
spin-3/2 case, models (J2J and have an exact SO (5) 
~ Sp(4) symmetry p|. This Sp(2AT) symmetry consider- 
ably simplifies the problem but may appear rather arti- 
ficial. However, we expect that, for generic and small in- 
teractions, the original SU(2) spin-rotational invariancc 
will be dynamically enlarged at sufficiently low energy 
0,0- A second reason to consider the Sp(2iV) symmet- 
ric model (J2J stems from the fact that the Sp(2iV) and 
SU(2) groups share the same center, the Z 2 group. More- 
over, the striking physical properties of the system rely 
on the existence of a Z^r symmetry which is also a sym- 
metry of the SU(2) model This %jq symmetry, that 
simply amounts to a global redefinition of the fermion 
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phase, is properly defined as the coset between the cen- 
ter I2N of the SU(2iV) group and the center Z2 of the 
Sp(2JV) or SU(2) one: Z N = Z 2N /Z 2 with 
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2N 
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the Z 2 symmetry being c a A — ► — c Q;! ;. The Zn symme- 
try, defined by Eq. J3Jl with n = 0, ...,7V— 1, provides an 
important physical ingredient not present in the F = 1/2 
case. The stabilization of a quasi-long-range BCS phase 
for F > 1/2 requires the spontaneous breaking of this 
Zjy symmetry since the singlet pairing 7r] is not invariant 
under this symmetry. Since Zn may be also viewed as 
a discrete subgroup of the global U(l) charge symmetry 
(see Eq. @), it is tempting to interpret, for F > 1/2, 
the breaking of Zn as the reminiscence of the sponta- 
neous global U(l) charge breaking that characterizes the 
BCS phase in higher dimensions. In contrast, if Zn is 
not broken, the BCS instability is suppressed and the 
leading superfluid instability, which has to be a singlet 
under the Zn symmetry, is a molecular object made of 
2N fermions. In the following, the delicate competition 
between these superfluid instabilities will be investigated 
by means of a low-energy approach. 

-Phase diagram. The low-energy effective field the- 
ory associated with Eq. J5J) is obtained, as usual, from 
the continuum description of the lattice electronic oper- 
ators in terms of right and left moving Dirac fermions: 
c a ,i/\/ao — * R a {x)e tkFX + L a {x)e~ lkpX , x = ia , ao be- 
ing the lattice spacing and kp is the Fermi momentum 
[7j. Away from half-filling (i.e. N atoms per site), it 
separates into two commuting density and spin pieces: 
H = H d + H s with [Hd,H s ] = 0. The U(l) density sec- 
tor is described by a bosonic field <& and its dual whose 
dynamics is governed by the free-boson Hamiltonian: 
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where v = v F [1 + (2V + UN{2N - l))/(iV7rw F )] 1/2 (v F = 
2tao sin(/c F ao) being the Fermi velocity) and K c i = 
[1 + (2V + UN(2N - l))/(N7rv F )r 1/2 are the Luttinger 
parameters. The conserved quantities in this U(l) sec- 
tor are the total particle number and current: Af 
J dx (RlR a + LlL a ) 
Jdx (~RlR a + LlL a ) 



y/2N/n J dx d x $ and J = 



y2iV /tt J dx d x Q respectively. 
For incommensurate fillings, the density degrees of free- 
dom are massless and display metallic properties in the 
Luttinger liquid universality class Q- All non-trivial 
physics is encoded in the spin part of the Hamiltonian: 
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where g\\ = -2{2V + NU)/N, g ± = 2(2V - NU)/N 
and we have neglected a velocity anisotropy. The Hamil- 
tonian JSJ) describes a SU(2iV)i conformal field theory 
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FIG. 1: Phase diagram of model Dashed lines denote 
cross-over lines whereas the solid; this results in strong con- 
straints on the possible phases line marks the phase transition 
in the Zjv universality class between phases II and III. 



(CFT) perturbed by a marginal current-current interac- 
tion. In Eq. JSJ, the currents Ir( L \, A = (1, ...,4/V 2 - 1), 
of the SU(2iV)i CFT have been decomposed into |j and 
_L parts I A = with respect to the Sp(2A) sym- 



R(L) 



metry of the lattice model (J2J. The currents Zj? 
l,..,N{2N + 1) generate the Sp(27V)i CFT symmetry 
and can be simply expressed in terms of the chiral 
Dirac fermions: K R = R^T^pRp, T a being the gen- 
erators of Sp(27V) in the fundamental representation. 
The remaining SU(27V)i currents are written as F L ±R = 
K^T^pRp, i = 1, ..,2N 2 - N - 1 (similar expressions hold 
for the left currents). The next step of the approach is to 
consider a description which singles out the Zn symme- 
try of the lattice model discussed above. To this end, 
we shall use the quantum equivalence 0: U(2N)i — > U(l) 
x Sp(2/V)i x Ztv, Zn being the parafermion CFT which 
describes self-dual critical points of two-dimensional Zn 
Ising models [IJ. This conformal embedding provides us 
with a non-perturbative basis to express any physical op- 
erator in terms of its density and spin degrees of free- 
dom which are described respectively by the U(l) and 
Sp(2A)i x Z N CFTs. The lattice Z N symmetry is then 
captured, within this low-energy approach, by an effec- 
tive 2D Zn model which is a generalization to TV > 2 of 
the standard Ising model. As in the N ~ 2 case, these 
Zn Ising models exhibit two phases described by order 
and disorder parameters and k = 1, .., N— 1, which 
are dual to each other by means of the Kramers- Wannier 
(KW) duality symmetry. This duality transformation 
maps the Zn symmetry, broken in the low-temperature 
phase ((crfc) ^ and (//&) = 0), onto a Zn symmetry 
which is broken in the high-temperature phase where 
((J>k) 7^ and ((Tfe) = 0. At the critical point, the theory 
is self-dual with a Zn X Zn symmetry and its univer- 
sal properties are captured by the Zn parafermion CFT 
0. In the simplest N = 2 case, there is a simple free- 
field representation of the unperturbed SU(4)i CFT in 
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terms of real fermions which has been extensively used 
in the context of two- leg ladders 0]. Introducing real 
fermions £ RL and £ R l,cl = 1, ...,5 to describe respec- 
tively the Z 2 and SO(5)i ~ Sp(4) x CFTs, the interacting 
part of Eq. © becomes Wj* = g^Ut) 2 +9^UUUl- 
The latter model has been studied recently to describe a 
SO(5) symmetric two-leg ladder [T^j. F° r generic N, the 
phase diagram of Eq. © can be elucidated by means of 
a two- loop renormalization group (RG) analysis. As de- 
picted in Fig. 1 it consists of three regions Region 
I is a generalized spin-density-wave (SDW) phase which 
is obtained when U and V are positive. In that case, 
ffj-jll ^ in the infrared (IR) limit and the interaction is 
marginal irrelevant. Up to a spin- velocity anisotropy, the 
low-energy properties of this phase are similar to that of 
the repulsive SU(2iV) Hubbard chain with 2N — 1 gapless 
spin excitations |l2l Il3j. In contrast, a spin gap opens 
in the two remaining regions which are distinguished by 
the Iin symmetry. In phase II, defined by U < and 
V > NU/2, the RG flow in the far IR limit is attracted 
along a special symmetric ray g\\ = g± = g* > where 
the interacting part of the Hamiltonian © can be rewrit- 
ten in a SU(2iV) invariant form: 

K nt = .9* (W + ArAl) = 9* lilt (6) 

The Hamiltonian ©, which governs the IR properties 
of phase II, takes the form of the SU(2iV) Gross-Neveu 
(GN) model which is an integrable massive field theory 
|l4j . It is instructive to have a simpler understanding 
of the spin-gap formation from the underlying Zjv Ising 
model. The Zjy and In symmetries that define the low- 
T and high-T phases of this Ising model admit a repre- 
sentation in terms of the fermions (R,L) a . Indeed, the 
Zn group © can be generated in this chiral fermion 
basis with help of the unitary operator U = c i7rA ^/ JV : 
UR{L)jA = e iw / N R(L) a . Similarly, the dual Z N sym- 
metry can be defined by U = e ™ J ' N : UR{L) a U^ = 
e^ '/ N R(L) a . The ground state of model © displays 
long-range order associated with the order parameter 
Tr( 5 ) = RlL^V^^. In phase II, we find that the 
Z2W = Z2 x Zyv symmetry remains unbroken while Z^v 
is spontaneously broken. The Zjy Ising model thus be- 
longs to its high-T phase and a spectral gap is formed. 
In the second spin-gapped phase (III) of Fig. 1, defined 
by V < and V < NU/2, the RG flow is now attracted 
along the asymptote: g\i = —g± = g* > 0. In that case, 
the interacting part of the IR Hamiltonian becomes 

K n * = 9* ( 7 |ffl 7 p, ~ ArAl) , (7) 

which can be recast as a SU(2iV) GN model © by 
means of a duality transformation T> on the fermions: 
VR^V- 1 = R(L) with R a = J a pR^ and L a = L a . 

This transformation acts on the currents as: -fffam = 
I\\R(L) an d Aril) = — (+)Ar(l) so tna ^ ^ m deed maps 



© onto ©. Besides the opening of a spectral gap, 
we thus find that phase IHjpossesses a hidden symme- 
try at low energy i.e. a SU(2N) symmetry generated 
by the dual currents {I^b,(l)^Ar(l))- ^ n ^ ac ^' one nas 
VUV- 1 = U so that V identifies to the KW duality of 
the Zjy Ising model. In phase III, the latter model thus 
belongs to its low-T phase and the In symmetry is spon- 
taneously broken whereas Z 2 jv = Z 2 x Zjy remains un- 
broken. In summary, the existence of these two distinct 
spin-gapped phases is a non-trivial consequence of higher- 
spin degeneracy and does not occur in the F = 1/2 case. 
The emergence of the spin-gap stems from the sponta- 
neous breakdown of the Z^r or In discrete symmetries. 
As we shall see now, these symmetries are central to the 
striking physical properties displayed by these phases. 

-Spin superfluidity. The low-energy properties of the 
spin sector of phase II can be extracted from the inte- 
grability of the SU(2iV) GN model ©. Its spectrum 
consists into 2N — 1 branches that transform in the 
SU(2iV) representations 0|. These eigenstates are la- 
belled by quantum numbers associated with the con- 
served quantities of the SU(2iV) low-energy symmetry: 
Qjf = Jdx {q L a = (1, N) and Q\ = 

J dx (P XL + P XR ), i = (1,...,N - I). Due to the 
Sp(2iV) symmetry of model (0, the Qjj numbers are 
conserved whereas the Q\ charges are only good quan- 
tum numbers at low energy. The spin spectrum in phase 
III can be obtained from the duality symmetry T> and 
consists into 2N — I branches which transform in the 
representations of the dual group SU(2N). The dual 
quantum numbers are now given by: Qjj = Qjj and 

= Jdx (P ±L + P ±R ) = Jdx (P ±L P ±R ) = J*. 
We thus observe that the low-lying excitations in phase 
III carry quantized spin currents in the "_L" direction. 
In this sense, the phase III might be viewed as a partially 
spin-superfluid phase. 

-Confinement. We shall now determine the nature 
of the dominant electronic instabilities of the different 
phases of Fig. 1. To this end let us consider oper- 
ators O raj which carry particle number n and current 
j- W(J),O n ,j] = n{j)O n .j. Using the Z 2 jv and Z 2JV 
generators, IA and Li, we find that O n j carry Z 2 ^r and 
Z 2 jv charges n and j respectively. In phase II, the full 
Z 2 jv = Z 2 x 1^ symmetry © is unbroken so that it 
costs a finite energy gap to excite states that either 
break the Z 2 or In symmetries. The dominant insta- 
bilities must thus be neutral under the Z 2 jv symmetry 
and the resulting order parameters O n ,j arc character- 
ized by 71 = mod 2N and j = mod 2. In particular, 
there is no quasi-long-range BCS order in phase II since 
the lattice singlet pairing operator irj carries a charge 
11 = 2 under the Z 2 n symmetry ©. The Z 2 n sym- 
metry thus confines the electronic charge to multiple of 
2Ne i.e. the leading supcrfluid instability in phase II 
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is a composite object made of 2N fermions. In this re- 
spect, the dominant order parameters in phase II are: 
P 2k r = LlR a and H^ = ^-^^...^L^.L^ 
which are, respectively, the 2k F component of the atomic 
density p and the uniform component of the lattice 
SU(2A)-singlet superconducting instability made of 2N 
fermions: U 2N Ui) = e ai - a2N c\ ..A ,. These orders 
are power-law fluctuating: (p\ k (a;)p 2 fc P (0)) ~ x~ Kd ^ N , 
and (ll 2 m (x)ll 2 N {0)) ~ x~ N /' K *. For K d < N, the 
leading instability is p2k F which gives rise to an atomic- 
density wave (ADW) phase whereas for Kd > N a 
SU(2iV) molecular-supcrfluid (MS) phase is stabilized 
(see Fig. 1) with order parameter IIq . The prop- 
erties of phase III are obtained from those of phase II 
with help of the duality symmetry V : (A/ <-> J ', Zjy <-> 
Z^,Kd <-> 1/Kd). Low-energy excitations in phase III 
carry now n = mod 2 and j = mod 2N since the 
symmetry 1*2 n = ^2 x ^at remains unbroken. We find 
now the confinement of atomic currents and the emer- 
gence of a quasi-long-range BCS pairing phase. Under 
the duality V symmetry, the ADW phase is mapped 
onto a BCS phase for K d > 1/N with order parameter 
= Ra^apL^p whereas the MS phase is mapped onto a 
molecular density- wave (MDW) phase with order param- 
eter n 2N — ,qi../3jv <j <j D D rt rt 

which emerges when K d < 1/N. The spontaneous 
breaking of the Ijq symmetry |j3J thus accounts for the 
emergence of the BCS superfluid phase and the spin- 
superfluidity phenomenon discussed above. The possible 
occurrence of two different superfluid phases II and III 
may be probed experimentally. Consider for example a 
gas trapped in an optical potential of length L, with an 
annular geometry and moving with tangential velocity 
V. This amounts to imposing a total particle current in 
the system J = 4NV/Vo where Vo = h/mL. In the su- 
perfluid phase III since the low energy excitations carry 
currents j = mod 2N we expect the total energy E(V) 
to display degenerate minima for quantized velocities : 
V n = nVo/2 irrespective of the value of the spin F. In 
contrast in the phase II where currents are unconfined, we 
expect the degenerate minima of E(V) at V n = nVo/2N. 

-The Zjy phase transition. The nature of the quantum 
phase transition between the two spin-gapped phases II 
and III can be determined through the duality symmetry 
V. On the self-dual line g ± = 0, i.e. 2V = NU, there is 
a separation of the Sp(2A) and Zjv degrees of freedom. 
Though the Sp(2A) sector remains gapfull when U < 0, 
the effective In Ising model is at its self-dual critical 
point and governs the phase transition. The Zjy quantum 
criticality for N = 2, 3 may be revealed by considering 
the ratio TZ N (x) = {^\{x)^))) N '' '/ 1 ^\x)Uf '(0)) 
which displays a power-law decay with a universal ex- 
ponent: TZ N (x) - a;-2iV(JV-i)/(^-i-2)_ For larger N the 

phase transition is non-universal. For N > 4, a strongly 
relevant perturbation is indeed generated which takes the 



form of the second thermal operator e? of the %jq CFT 
with scaling dimension 12/(A + 2) [9|. The resulting 
model is intcgrable and the transition is either of first- 
order or in the U(l) universality class depending on N 
and the sign of the coupling constant of e 2 0. In the 
N = 4 case, i.e. a special case of the Ashkin- Teller model, 
the Z4 criticality can emerge with the introduction of the 
interaction AJ^Pqq ^-Poo, i) 2 which can eliminate the op- 
erator 62 by a fine-tuning of A. 

-Mott phases. At the commensurate 1/2N filling, i.e. 
one atom per site, an umklapp term ~ cos (V87rAT$) 
is generated in the density sector and becomes relevant 
when Kd < 1/N leading to the opening of a density 
gap. We further distinguish between three different Mott 
phases 0| . The first one lies in the SDW region I and is 
qualitatively similar to the one encountered in the pure 
SU(2A) Hubbard chain 0. The two others have a spin 
gap and can be distinguished with respect to the confine- 
ment properties of atomic currents. In the region II we 
find for large enough V a 2k F -ordered ADW and spin- 
Peierls ordering with a 2N ground-state degeneracy. In 
the region III the cross-over line between BCS and MDW 
phases identifies to the Mott transition line. At filling 
1/2N, while the BCS phase remains, the MDW regime 
locks into a 2Nk F MDW and displays a dimerized bond 
ordering for all N. This Mott phase is unusual since there 
is no one-particle density long-range fluctuation due to 
the confinement of atomic currents. Remarkably enough 
we find that this MDW Mott phase is the only gapped 
phase directly connected to the BCS superfluid phase. 
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